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High-energy quark-quark scattering and the eikonal approximation 
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The high-energy quark-quark scattering amplitude is calculated first in the case of scalar QCD, using Fradkin's 
approach to derive the scalar quark propagator in an external gluon field and computing it in the eikonal approx- 
imation. The results are then extended to the case of "real" (i.e., fermion) QCD. The high-energy quark-quark 
scattering amplitude turns out to be described by the expectation value of two lightlike Wilson lines, running 
along the classical trajectories of the two colliding particles. Interesting analytic properties of the high-energy 
quark-quark scattering amplitude can be derived, going from Minkowskian to Euclidean theory: they could open 
the possibility of evaluating the high-energy scattering amplitude directly on the lattice. 
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It is known that for a class of soft high-energy 
scattering processes, i.e., elastic scattering pro- 
cesses at high squared energies s in the center 
of mass and small squared transferred momen- 
tum t (that is s (X) and |t| <C s, let us say 
\t\ < 1 GeV^), QCD perturbation theory cannot 
be safely applied, since t is too small. Elabo- 
rate procedures for summing perturbative contri- 
butions have been developed Q, even if the 
results are not able to explain the most relevant 
phenomena. 

A non-perturbative analysis, based on QCD, of 
these high-energy scattering processes was per- 
formed by Nachtmann in He studied the 
s dependence of the quark-quark (and quark- 
antiquark) scattering amplitude by analytical 
means, using a functional integral approach and 
an eikonal approximation to the solution of the 
Dirac equation in the presence of an external non- 
Abelian gauge field. 

In a previous paper we proposed an ap- 
proach to high-energy quark-quark (and quark- 
antiquark) scattering, based on a first-quantized 
path-integral description of quantum field the- 
ory developed by Fradkin in the early 1960's 
In this approach one obtains convenient expres- 
sions for the full and truncated-connected scalar 
propagators in an external (gravitational, elec- 
tromagnetic, etc.) field and the eikonal approx- 
imation can be easily recovered in the relevant 
limit. Knowing the truncated-connected propa- 



gators, one can then extract, in the manner of 
Lehmann, Symanzik, and Zimmermann (LSZ), 
the scattering matrix elements in the framework 
of a functional integral approach. This method 
was originally adopted in 1^ in order to study the 
four-dimensional Planckian-energy gravitational 
scattering. 

One starts, for simplicity, with the case of 
scalar QCD, i.e., the case of a spin-0 quark 
(described by the scalar field 0) coupled to a 
non-Abehan gauge field A^" = Ai^Ta, Ta (a = 
1, . . . , N"^ — 1) being the generators of the Lie al- 
gebra of the colour group SU{Nc). 

In the Fradkin's approach Q one writes the 
propagator for the scalar field as a functional 
integral of the first quantized theory |^ . 

The Green function (Feynman propagator) for 
the scalar field (f) in the external metric g^i, and in 
the external non-Abelian gauge field admits 
a representation in terms of a functional integral 
over trajectories H M: 



G{y,x\g,A) = 



dv 
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where 
L(X,X) = 
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The matrix 1 is the Nc x Nc unity matrix in 
the colour space and "P" means ^^path-ordered'' 
along the "histories" X^(r). Moreover g{X) = 
det(5^^(X)), so that g{X) = — 1 in a flat space- 
time. The propagator (1) is already connected, 
i.e., vacuum diagrams have been already divided 
out. Yet, in order to derive the scattering matrix 
elements following the LSZ approach, we need 
to know the on-shell truncated-connected Green 
functions, which are obtained from the connected 
Green functions by removing the external legs cal- 
culated on-shell. In Ref. Q it is shown how 
to compute the truncated-connected propagator 
Gt{p,p'\A), starting from Eq. (1), in the so- 
called eikonal approximation, which is valid in the 
case of scattering particles with very high energy 
(i? = p° ~ IpI ^ m) and small transferred mo- 
mentum q = p' — p (i.e., \t\ <C E, where t = q^). 
In this limit the functional integral can be eval- 
uated by means of a saddle-point approximation 
in which the classical trajectory is computed to 
the lowest non-trivial order, that is a straight 
line (in the Minkowski space-time) described by 
a classical particle (of mass m) moving with four- 
momentum pf^ : 



m 



(3) 



One thus finds, after properly generalizing the 
results of to the case of an external non- 
Abelian gauge field, the following expression 
for the truncated-connected propagator in the 
eikonal approximation: 



Gt{p,p'\A)^2E / d^be 



Pexp 
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1 (4) 



As said before, q = p' — p is the transferred mo- 
mentum. In one must not include the com- 
ponent of which is parallel to pf^. For ex- 



ample, if p^ 



P 



{E,E, 0,0), one has that 



p^ ~ p'^ ~ 2E and p_ ~ p'_ ~ 0, where the gen- 
eral notation A+ = A" + A\ yl_ = A° - A^ 
has been used for a given four-vector A''. In 
this case one must take: d^b — d^htdb^, where 
ht = {by, bz) is the component of b^ in the trans- 
verse plane (y, z) (while the two light-cone coor- 



dinates, x+ = t+x and x- = t — x, are sometimes 
called longitudinal coordinates) . 

Knowing the truncated-connected propaga- 
tors, one can then extract, using the reduction 
equations of LSZ in the framework of a func- 
tional integral approach, the scattering amplitude 
of two high-energy spin-0 quarks, in the limit 
s ^ oo and t ^ s. The results can then be easily 
extended to the case of physical interest, i.e., the 
case of "real" (fermion) QCD. 

The quark-quark scattering amplitude, at high 
squared energies s in the center of mass and small 
squared transferred momentum t (that is s — > oo 
and 1^1 <^ s, let us say \t\ < 1 GeV^), turns out 
to be described by the expectation value of two 
lightlike Wilson lines, running along the classical 
trajectories of the two colliding particles [|j 

In the center-of-mass reference sys- 
tem (c.m.s.), taking for example the initial tra- 
jectories of the two quarks along the x^-axis, the 
scattering amplitude has the following form [ex- 
plicitly indicating the color indices («, j, . . .) and 
the spin indices (a, /3, . . .) of the quarks] 

Mf^ = {lp^a{p'l)ll^k'({P2)\M\lpjl3{pi)'4'ls{P2)) 



■ Saf36.ys ■ 2s / d^zte''^-^' X 



x{[Wiizt) - ll,j[W2{0) - l]ki) , 



(5) 



where q = (0,0, q), with t ^ q^ = — q^, is the 
tranferred four-momentum and zt = (0, 0,zt), 
with Zt ~ (z^,z"^), is the distance between the 
two trajectories in the transverse plane. The ex- 
pectation value if (A)) is the average of f{A) in 
the sense of the functional integration over the 
gluon field A'^ (including also the determinant of 
the fermion matrix, i.e., det[i7^£)^ — m], where 
_ gn igAf^ is the covariant derivative) ^ 
g. The two lightlike Wilson lines Wi{zt) and 
14^2(0) in Eq. (5) are defined as 



Wi{zt) = Pex^p 



W2{0) = Pexp 



/+00 
A^,{zt + PlT)p{dT 
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where P stands for ^^path ordering" and A^^ = 
A^T'^; pi ~ {E,E,Qt) and p2 =^ {E,-E,Qt) 
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Fig. 1. The space-time configuration of the two 
lightlike Wilson lines Wi and W2 entering in the 
expression (5) for the high-energy quark-quark 
elastic scattering amplitude. 



are the initial four-momenta of the two quarks. 
The space-time configuration of these two Wilson 
lines is shown in Fig. 1. 

Finally, in Eq. (5) is the fermion-field 
renormalization constant, which can be written 
in the eikonal approximation as 



{TT[w,m) 



N, 



■{Tr[W2m) 



(7) 



In a perfectly analogous way one can also derive 
the high-energy scattering amplitude in the case 
of the Abelian group C/(l) (QED). The resulting 
amplitude is equal to Eq. (5), with the only obvi- 
ous difference being that now the Wilson lines Wi 
and W2 are functionals of the Abelian field A^^ 
(so they are not matrices). Thanks to the sim- 
ple form of the Abelian theory (in particular to 
the absence of self-interactions among the vector 
fields), it turns out that it is possible to explicitly 
evaluate (at least in the quenched approximation) 
the expectation value of the two Wilson lines [Q : 
one finally recovers the well-known result for the 
eikonal amplitude of the high-energy scattering 
in QED [| § 0. 

In what follows, we shall deal with the quantity 



9M(i3,kl){s; t) 
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x{[WM)-'^h[W2{Q)-l]ki) , (8) 

in terms of which the scattering amplitude can be 
written as 

Mf, = {ll},a{p'l)^k-i{p'2)\M\ll)jl3{pi)^u{p2)) 

~i-2s- 5ap&^s ■ gM{ij,ki) t) ■ (9) 

The quantity gM{ij,ki) ^) depends not only on 
t = — q^, but also on s. In fact, as was pointed out 
by Verlinde and Verlinde in it is a singular 
limit to take the Wilson lines in (8) exactly light- 
like. A way to regularize this sort of "infrared" 
divergence (so called because it essentially comes 
from the limit m — > 0, where m is the quark 
mass) consists in letting each line have a small 
timelike component, so that they coincide with 
the classical trajectories for quarks with a finite 
mass m (see also Ref. [|l2| for a discussion about 
this point). In other words, one first evaluates the 
quantity gM{ij,ki) iP'^ ^) '^'-"^ Wilson lines along 
the trajectories of two quarks (with mass m) mov- 
ing with velocity (3 and —(3 (0 < /3 < 1) in the 
x^-direction. This is equivalent to consider two 
infinite Wilson lines forming a certain (finite) hy- 
perbolic angle x ii^ Minkowski space-time. Then, 
to obtain the correct high-energy scattering am- 
plitude, one has to perform the limit /3 — > 1, that 
is X ^ 00, in the expression for gM(ij,fci) (/?; t)'- 

Mf, = {'ll;,a{p'l)^k-i{p'2)\M\lpjfi{pi)^l5{p2)) 

^ -i-2s- SapS^s ■ gMUj,ki)iP 1' (10) 

s—^ao 

Proceeding in this way one obtains a In s depen- 
dence of the amplitude, as expected from ordi- 
nary perturbation theory and as confirmed by the 
experiments on hadron-hadron scattering pro- 
cesses | [|. 

For deriving the dependence on s one exploits 
the fact that both (3 and x ^-re dependent on s. 
In fact, from E = m/ \J\ — (fi and from s = 4£^^, 
one immediately finds that 



/3 = \ 1 



4m^ 



(11) 



By inverting this equation and using the relation 
/3 — tanh?/;, where '0 is the hyperbolic angle [in 
the plane (a;°,a;^)] of the trajectory of Wi, we 
derive that 



s = 4m^ cosh^ = 2m^(coshx + 1) 



(12) 
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Therefore, in the high-energy hmit s — > oo (or 
/3 — > 1), the hyperbohc angle x = SV' is essentially 
equal to the logarithm of s (for a finite non-zero 
quark mass m): 



In s 



(13) 



In Sect. 3 of Ref. |T|] we have followed this pro- 
cedure to explicitly evaluate the second member 
of (10) up to the fourth order in the expansion 
in the renormalized coupling constant gj^: the re- 
sults so derived are in agreement with those ob- 
tained by applying ordinary perturbation theory 
to evaluate the scattering amplitude up to the 
order 0(g|) § §. ^ 

The direct evaluation of the expectation value 
(8) is a highly non-trivial matter and it is 
strictly connected with the ultraviolet proper- 
ties of Wilson-line operators ||l^. Some non- 
perturbative approaches for the calculation of (8) 
have been proposed in Refs. [Q and |l^]. At 
the moment, the only non-perturbative numer- 
ical estimate of (8), which can be found in the 
literature, is that of Ref. [|l5| (where it has been 
generalized to the case of hadron-hadron scatter- 
ing): it has been obtained in the framework of 
the stochastic vacuum model (SVM). 

In some recent papers jlj] we have pro- 
posed a new approach, which consists in adapt- 
ing the scattering amplitude to the Euclidean 
world: this approach could open the way for the 
direct evaluation of the scattering amplitude on 
the lattice. More explicitly, we have shown that 
the expectation value of two infinite Wilson lines, 
forming a certain hyperbolic angle in Minkowski 
space-time, and the expectation value of two in- 
finite Euclidean Wilson lines, forming a certain 
angle in Euclidean four-space, are connected by 
an analytic continuation in the angular variables. 

Let us consider the following quantity, defined 
in Minkowski space-time: 



9m{pi,P2; t) = ^ [ d^zte*i-^'x 
x{[W^{zt)-lUW2{0)-l]ki) 



(14) 



where pi and p2 are the four-momenta [lying (for 
example) in the plane (a:°,a;^)], which define the 
trajectories of the two Wilson lines Wi and W2 



{Af^ — A'^T'^ and m is the fermion mass^ 
Wi{zt) 



Pexp 



11^2(0) = Pexp 
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A,i?lrfldr 



m 



Z\fi/ in Eq. (14) is defined as 



^(Tr[I^i(0)]) = ^(Tr[I^2(0)]) 



(15) 



(16) 



This is a sort of Wilson-line's renormalization 
constant: as shown in Ref. j^], Zy^ coincides with 
the fermion renormalization constant Z^ in the 
eikonal approximation. 

By virtue of the Lorentz symmetry, we can de- 
fine pi and p2 in the c.m.s. of the two parti- 
cles, moving with speed (5 and —(5 along the x^- 
direction: the hyperbolic angle x [iii the plane 
(a;°,a;^)] between the two trajectories of W\ and 
Wi is thus given by x = 2 arctanh/?. 

In an analogous way, we can consider the fol- 
lowing quantity, defined in Euclidean space-time: 

gE{piE,P2E] t) = I d^zte'*!""* X 

y{[WiE{ztE) - 1]^J[W2E{0) - l]kl)E , (17) 

where zts = (21,22,^3,^4) = (0,Zt,0) and qe = 
(0,q, 0) (so that: q% — = -t). The expec- 
tation value {.■.)e must be intended now as a 
functional integration with respect to the gauge 



variable A 



(E) 



A^E)aj,a ^Yie Euclidean theory. 



The Euclidean four-vectors piE and p2E [lying 
(for example) in the plane (xi, 2:4)] define the tra- 
jectories of the two Euclidean Wilson lines Wie 
andV(^2_B- .^wfiinEq. (17) is defined analogously 
to Zw in Eq. (16): 

ZwE ^ ^(Tr[Tyis(0)]) = ^(Tr[I^2£(0)]) .(18) 

Denoting by the Euclidean angle in the plane 
(xi,a;4) between the trajectories of the two Eu- 
clidean Wilson lines Wie and W2E^ one finally 
finds the following relation between the ampli- 
tudes 5m(x; t) and .g£;(6'; t) 



gAiix; t) — '9M{i6; t) 



9e{0; t) 



9E{e- t) ■ 

9E{-ix; t) = 3m(x; t) 



(19) 



5 



This relation of analytic continuation has been 
proven in Ref. for an Abelian gauge theory 
(QED) in the so-called quenched approximation 
and for a non-Abelian gauge theory (QCD) up 
to the fourth order in the renormalized coupling 
constant in perturbation theory. In Ref. we 
have generalized the results of Ref. giving the 
rigorous proof of the above-mentioned relation 
of analytic continuation for a non-Abelian gauge 
theory with gauge group SU{Nc) [as well as for 
an Abelian gauge theory (QED)]. The approach 
adopted in Ref. consisted in explicitly eval- 
uating the amplitudes gM and 5b in some given 
approximation (such as the quenched approxima- 
tion) or up to some order in perturbation theory 
and in finally comparing the two expressions so 
obtained. Instead, in Ref. [|l6| we have given a 
general proof, which essentially exploits the rela- 
tion between the gluonic Green functions in the 
two theories. 

Therefore, it is possible to reconstruct the 
high-energy scattering amplitude by evaluating 
a correlation of two infinite Wilson lines forming 
a certain angle 9 in Euclidean four-space, then 
by continuing this quantity in the angular vari- 
able, 9 —>■ ~ix, where x is the hyperbolic angle 
between the two Wilson lines in Minkowski space- 
time, and finally by performing the limit x ^ oo 
(i.e., f} — > 1). In fact, the high-energy scattering 
amplitude is given by 

Mfi = {'^|)^a{p'l)i^k-1{J>'2)\M\■^l)jf^{pl)^ls{p2)) 

- -i-2s- SapS^s ■ gnix oo; t) . (20) 

s — >oo 

The quantity gniXi t)i defined by Eq. (14) in the 
Minkowski world, is linked to the corresponding 
quantity gE{9; t), defined by Eq. (17) in the Eu- 
clidean world, by the analytic continuation (19) 
in the angular variables. The important thing to 
note here is that the quantity gE{9] t), defined in 
the Euclidean world, may be computed non per- 
turbatively by well-known and well-established 
techniques, for example by means of the formu- 
lation of the theory on the lattice or using the 
SVM |l|. In all cases, once one has obtained 
the quantity gE{9] t), one still has to perform 
an analytic continuation in the angular variable 
9 — > ~ix, and finally one has to extrapolate to 



the limit x ^ ^ (i-G-, — * !)• For deriving the 
dependence on s one exploits the fact that both 
P and X are dependent on s, according to Eqs. 
(11), (12) and (13). 

Of course, the most interesting results are ex- 
pected from an exact non perturbative approach, 
for example by directly computing (6*; t) on the 
lattice: a considerable progress could be achieved 
along this direction in the near future. 
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DISCUSSIONS 

L.N. Lipatov, INP, St. Petersburg (Russia) 
What is the relation of your work with the BFKL 
results (LLA for high-energy scattering)? In 
1986 the quark quark scattering amplitudes with 
the elastic and quasi-elastic unitarity were con- 
structed. Are they related with your eikonal am- 
plitude? 

E. Meggiolaro 

It is known that "soft" hadronic collisions at high 
energy can be described by the Pomeron exchange 
model. However, a more basic understanding of 
Pomeron exchange in terms of the underlying 
quark-gluon dynamics (QCD) remains an out- 
standing open problem. Most theorists now agree 
that the properties of the Pomeron can be some- 
how derived in QCD from multi-gluon exchange 
and the physical picture commonly accepted is the 
one where the Pomeron couples to single partons: 
yet, it is still unclear how to develop a reliable cal- 
culational framework based on this idea. The first 
calculations of the near forward parton parton 
scattering amplitudes were done to the lowest or- 
ders of perturbation theory (PT) in the leading 
logarithmic approxim,a,tion (LLA), or using the 
elastic and quasi elastic unitarity. In particular, 
the high-energy asymptotics of the scattering am- 
plitude was found in the LLA, using the results 
of one-loop calculations in order to resum the 
leading logarithmic corrections to all oders of PT. 
However, in spite of the considerable progress so 
achieved, there was no regular way to resum and 
control non-leading logarithmic corrections to the 
scattering amplitude. 

The approach described in my work enables us 
to take into account systematically these non- 
leading corrections. It starts from the first 
principles of QCD and derives the high-energy 
quark-quark scattering amplitude in a fully non- 
perturbative way, using a functional integral ap- 
proach and the so-called eikonal approximation 
to write the truncaied -connected quark propa- 
gator in an external gluon field. A compar- 
ison of this approach to the standard PT re- 
sults has been performed up to the one-loop level 
(fourth order in the renormalized coupling con- 



stant): the two approaches give the same results. 
An explicit comparison with the BFKL Pomeron 
model is not yet available. It would be an ex- 
tremely important result to derive the features 
of the Pomeron model starting from my non- 
perturbative approach, based on first principles of 
QCD: this is matter for future work. 



